UCD-97-16 



hep-ph/9708264 



MESSENGER SNEUTRINOS AS COLD DARK MATTER 

Tao Han and Ralf Hempfling 

Davis Institute for High Energy Physics, Department of Physics 
University of California, Davis, CA 95616 USA 

a^ 

< 

^ ■ Abstract 

\^ . In models where supersymmetry breaking is communicated into the visible sector 

via gauge interactions the lightest supersymmetric particle is typically the gravitino 
■ which is too light to account for cold dark matter. We point out that the lightest 

, messenger sneutrinos with mass in the range of one to three TeV may serve as cold 

0^ ' dark matter over most of the parameter space due to one-loop electroweak radiative 

corrections. However, in the minimal model this mass range has been excluded by the 
Oh, direct dark matter searches. We propose a solution to this problem by introducing 

^ ^ terms that explicitly violate the messenger number. This results in low detection rate 

1^ I for both direct and indirect searches and allows messenger sneutrinos to be a valid 

"T^ \ dark matter candidate in a wide region of SUSY parameter space. 
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1 Introduction 



Theories in which the effects of super symmetry (SUSY) breaking are introduced into the "visible" 
sector via the standard model (SM) gauge interactions have recently received considerable atten- 
tion ||2|,^^,^,^^ . One of the most attractive features of such models is the natural explanation for the 
smallness of the SUSY-contributions to flavor-changing neutral current (FCNC) phenomena both in 
the quark and lepton sector as a result of a strongly constrained sparticle spectrum Q. A complete 
model of gauge mediated supersymmetry breaking requires three sectors: the "visible" sector con- 
taining the minimal supersymmetric standard model (MSSM); the "secluded" sector^ responsible for 
SUSY breaking and the messenger sector responsible for the communication of SUSY breaking effects 
into the visible sector.^ For the purpose of low energy phenomenology there is no need to specify 
any details about the secluded sector which in the minimal version provides only a single parameter 
(F) as explained below. The messenger sector is phenomenologically more important since it provides 
all soft SUSY breaking parameters for the MSSM. One generally assumes, that the messenger sector 
contains only complete SU(5) multiplets ( i.e. pairs of 5s and 10s) in order to naturally maintain 
gauge coupling unification (r-bottom unification on the other hand is problematic). The absence 
of a Landau-pole below the GUT scale limits the number of additional representations such that 

= A^5 + 3NiQ < 5 (this value depends somewhat on the messenger masses assumed to be universal 
and denoted by M). The general form of the messenger superpotential is given by 

VF = ^A$5$$, (1.1) 

where $ C 5, 10 denotes the messenger superfields, S is an MSSM singlet belonging to the secluded 
sector. In order to suppress FCNC effects we have to forbid any renormalizable interactions between 
messenger fields and MSSM matter fields by means of a symmetry. Hence, all the soft SUSY breaking 
MSSM parameters are determined by only two free parameters: M = X^{S) and A = F/{S) with 
F = (dWg/dS) (where Ws is the superpotential of the secluded sector), and a discrete choice for 
sign(/i) and A^. From here all the low energy parameters are obtained by renormalization group 
evolution In addition, there is the Higgs mass parameter |//| determined in the standard fashion 
by imposing radiative electroweak symmetry breaking. The /i parameter is not a soft SUSY breaking 
parameter and a mechanism has to be introduced in order to give it a phenomenologically allowed 
value Q. Such a mechanism typically weakens the predictability of the model by allowing the soft 
SUSY breaking Higgs mixing parameter Bq, or the ratio of the Higgs vacuum expectation values (vev) 
tan (3, to become a free parameter. 

One of the most important features of gauge mediated models is the fact that the gravitino is very 
light (its mass decreases with the square of SUSY breaking scale) while its coupling to the SUSY 



particle spectrum is enhanced ||10|. Consequently, the gravitino as the lightest SUSY particle can be 



an interesting warm [10| and mixed p] DM candidate, and the lightest neutralino, often heavier than 



the gravitino, is unstable and no longer a cold dark matter (CDM) candidate [pT|p!^ ]. The goal of 

"This sector replaces the so-called "hidden" sector in conventional supergravity (SUGRA) models. The fundamental 
difference from SUGRA models is the existence of renormalizable interaction in addition to gravity and, hence, a much 
larger ratio of the visible sector SUSY breaking scale over the secluded sector SUSY breaking scale. 

''Recently, there have also been attempts to combine messenger sector and secluded sector 
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this paper is to look for an alternative CDM candidate in the gauge mediated models. In section 2, 
we first reanalyze the properties of a messenger sneutrino as a CDM candidate. We then discuss the 
feasibility for its direct and indirect detections in section 3. Our conclusions are presented in section 4. 
Calculations and formula are relegated to three appendices. 



2 Messenger Sneutrinos as Cold Dark Matter 

Without an interaction of MSSM and messenger sector in the superpotential, the lightest messenger 
particle is stable and a possible CDM candidate. Several such candidates were investigated in Ref. 
In this paper, we will focus on the electrically neutral components of chiral supermultiplets, 5 and 5 
which was found to be the most promising possibility Q. It carries the same gauge quantum numbers 
as the MSSM scalar neutrinos and will be referred as messenger sneutrinos. We denote the electroweak 
doublets contained in 5 and 5 by $ and ^, respectively, and define ^± = (<I> it ^^)/^/2. In this basis 
the D-terms for <I> and <!• are 

= ($|^_T°$_ + $lr''$+) , (2.1) 
D' = |- + $L^>+) , (2.2) 
and the mass matrix can be written as 

v-i^i^i)( ™,f,+^.,„, «^"<rrs<°'>) m ■ M 



where M is the messenger mass scale. We have assumed that F is real and the electroweak indices 
have been suppressed. Without loss of generality we set F > 0. 

After electroweak symmetry breaking the D terms acquire a non-zero vev that lifts the mass de- 
generacy of the neutral and charged components of <I>±. The tree- level masses of the lighter mass 
eigenstates (p^ {Q = 0, —1 is the electric charge index) are M 



m^Q = M^- + (r3 - Q sin2 ew)^rni cos^ 2/3 . (2.4) 

In the SUSY limit (i.e. F = 0), the mass squared splitting of the neutral {Q = 0) and the charged 
[Q = —1) components is of order the Z mass squared m^. With SUSY-breaking and in the limit 
F » m^, we see that 

4 

m^- - m^o = sin2 29w cos^ 2/3 > . (2.5) 

This result looks promising: the lightest messenger particle is neutral. Furthermore, in a large part 
of the interesting parameter space the decay — > (j)^ff' is sufficiently fast in order not to disturb 
big bang nucleosynthesis (BBN). Unfortunately, as pointed out in Ref. 0, this suppression does 
in general not apply to the one-loop radiative corrections due to the SUSY breaking F-term. These 
corrections were found to be dominant and have the wrong sign over most of the parameter space, 
rendering the lightest messenger sneutrino (j)^ as a CDM candidate questionable. However, we find 
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that there are large one-loop corrections due to electroweak interactions even in the limit of unbroken 
SUSY. Numerically, we find that to a good approximation (obtained in the limit M ^ niz, F = and 
tan /? = 1), 

m^- — m^o ~ ^"^z — GeV . (2.6) 



This result dominates over the tree-level term of Eq. (2.5) and is only slightly modified by one-loop 



SUSY breaking effects for F ^ 0. Consequently, the decay occurs long before BBN and poses no 
significant constraint on the parameter space. The detailed results for the various one-loop diagram 
calculations are given in Appendix A. 

We now turn to the CDM relic density obtained from a freeze-out calculation of the messenger 
sneutrinos. The relic density of any particle is governed by the thermal average of the mean annihilation 
rate times the relative velocity, {av) at the freeze-out temperature, Tf. Typically, one finds that 
Tj ~ m^o /20 ^ THfj,- — m^o which implies that both charged and neutral components are present in 
the thermal bath at Tf. Hence, the freeze-out calculation should be most properly performed in the 
limit of symmetric phase for the four components {(jP , (j^'^)- We also ignore the gauge boson masses, 
the error in this approximation is of order vn^jM'^ and vanishes in the case that Ty is larger than 



the critical temperature of the electroweak phase transition. We follow Ref. |13| for the freeze-out 
calculation and our results are summarized in Appendix B. 

In Figure 1, we present our results for the relic density (0/i^) versus m^o for r = M/m^o = 1 
(solid curve) and r = 10 (dashed curve). The relic density for r = 1 is increased by about 50% for 
r = 10, and would go up by a factor of 2 for r 00. For comparison, we have also plotted the results 
calculated in the broken phase by considering only the lightest messenger sneutrino annihilation. 
The rather large discrepancy can be understood as follows: the broken phase calculation essentially 
assumes that all annihilation cross sections are of the same order, characterized by (lP'^(fP annihilation. 
However, notice that a{(f)4> anything) = (t((/)^(/)^ — > anything) = 0; and a"(</)^^(/>^ anything) is 
suppressed by one power of a; = Tf/m^i). As a result, the relic density calculated in the broken phase 
is underestimated by about a factor 4. Generally speaking, as see in Fig. 1, the sneutrino component 
of the SU(2) doublet with a mass in the TeV range could be a good CDM candidate. Imposing that 
the universe is not overdosed by the CDM particles ( i.e. Jl/i^ < 1) yields the upper bound on the 
mass of the lightest messenger sneutrino 

m^oOTeV. (2.7) 



For m^o < 1 TeV it seems inadequate to neglect Higgs and the squark masses. However, from Eq. ( |B.5D 
we see that the Higgs and matter fields are only produced in s-channel annihilation via intermediate 
gauge bosons. Hence, their rates arise at first order in x and only contribute a few percent. 



3 CDM Particle Detection 

The most direct way of detecting the existence of dark matter particles would be to observe the DM 



particle- nucleus scattering by recording the nuclear recoil in a detector [12|. It is known [14| that the 
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Figure 1: The relic density of <^ as a function of m^o in the limit that the mass of all MSSM particles 
are much lighter than m^o. The two upper curves are for our symmetric phase calculation and the 
lower ones are for the broken phase. 



cross section for the sneutrino-nucleus scattering via the Z exchange is 



a 



2tt {m^o + ruA 



[A-2{l-2Wew)Zf, (3.1) 



where Gp is the /i decay constant, is the nuclei mass, A the atomic mass number and Z the 
atomic number. The current direct searches put an upper bound on the cross section. In fact, a scalar 



DM particle with mass 3 TeV or less has been ruled out at a 90% confidence level [15| assuming it 
accounts for about 30% or more of the local galactic halo density 0.3 GeV/cm^. This would lead to 
the exclusion of the whole interesting region of the m^o parameter space in Fig. |l] and results in an 
unsatisfactory solution for the CDM issue in the minimal version of the gauge mediated model. 

To provide a solution to this problem, we propose to introduce a tree-level interaction of the MSSM 
Higgs sector and messenger sector. Consider the following superpotential (W) and soft SUSY breaking 
terms {V) 

W = \N^Hu + ]^mNN^ , 

V = \AxNmu + \Bxm^N^ + (3.2) 

where is a gauge singlet, the Higgs doublet coupled to up- type quarks and A\ and B\ the soft 
SUSY breaking parameters. The inclusion of an additional singlet N is motivated by the attempt to 
understand the problem of the ;U-parameter Q, namely, why the SUSY invariant /i-parameter is of 
the same order as the SUSY breaking scale in the visible sector. Eq. ( |3.2| ) may also have an effect on 
electroweak symmetry breaking |l^ which deserves further scrutiny. If both A and m^v are non-zero, 
then the messenger number is explicitly broken. However, the model still respects a Z2 parity which 
guarantees the stability of the lightest messenger particle. 
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After EW symmetry breaking, N is no longer a mass eigenstate and mixes with the neutral compo- 
nent of If we assume that rrij^ > M and )^{Hu)'^ <^ mf^ — M^, then the mixing is small so that our 
analysis in Sec. 2 is still valid. On the other hand, this mixing has another phenomenologically im- 
portant consequence: it leads to an additional mass splitting of the CP-even and CP-odd components 
of the neutral fields 



m 



mcp 

—even 



N 



2mj,o 



Bxm 



m 



N 



N 



M 



2Ay 



[m 



N 



N 



(3.3) 



BBN constraints are satisfied if this mass difference is larger than a few MeV. For instance, for 



M,m^o = 1 TeV, = 3 TeV and Bx,Ax = 100 GeV, Eq. implies A > 0(0.1). Since the 

DM particle coupling to the Z boson requires a CP-even and CP-odd transition, such a large mass 
difference would prevent DM particle-nuclei scattering via single Z exchange from happening. This 
is because the initial kinetic energy in the c. m. frame for DM particle-nuclei scattering is typically 
of order 0.1 MeV (assuming the DM particle velocity is about 10~^c) and is much smaller than the 
CP-even and CP-odd mass difference. The dominant contribution to DM particle-nuclei scattering 
would therefore be the spin-independent Higgs exchange through quarks and quark loops. In analogy 
to the neutralino-nucleon scattering [^] , we find that the cross section for sneutrino-nucleus scattering 
via Higgs exchanges is 



'^nucl 



4tt {m^o + 



[fpZ + fniA-Z)f, 



(3.4) 



where fp ~ are the effective sneutrino-nuclei coupling via Higgs exchanges. They are given in 
Apendix C. Following the discussion in Ref. [^], we estimate the direct search event rate {Rd) for 
scalar neutrino-nuclei scattering by 



Rd 



1.8 X lO^^GeV^ C f^nuci 



(3.5) 



kg day m^ouiA 
where C describes the suppression due to the nuclear form factor |15]. 
Another way of detecting the DM particle is by the large scale neutrino detectors via annihilation 



products from DM particles inside the sun/earth. Following the arguments in Ref. [19|, the expected 
event rate in a neutrino detector induced by z/^ from the (fp'^cjp annihilation in the sun may be estimated 
by 



2.65 X 1039 m^o 
10^ m^ yr CeV 



Cnucl 



cm^ 



(3.6) 



where is the relative abundance of element A with respect to hydrogen in the sun. With the 
nuclear form factor ( included, we find that the leading contributions to the DM particle capture are 
from the elements iron, oxygen and helium, rj in Eq. (|3.6| ) is the neutrino escape probability in the 



sun parameterized by |1£] 



1 + 



m. 



5.2 TeV 



(3.7) 
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Figure 2: The event rate for DM particles consisting of messenger sneutrinos in direct and indirect 
searches as a function of their mass, niM. 



which significantly suppresses the ly^ flux for m^o > 1 TeV. The background rate from the flux of 
atmospheric neutrinos coming from a pixel around the sun is estimated to be ]l9[ 



^ 10^ m2 yr VTevj 



The scattering cross section given by Eq. ( |3.4| ) is rather small. For a 1 TeV scalar neutrino scattering 
off a nucleon, it is typically of order 10~^^ ~ 10"^"^ cm^, depending on SUSY parameters. Comparing 
the direct and indirect detection rates given in Eqs. (3.5) and (3.6), we see that the indirect search 



seems to be more suitable for a heavier DM particle. Especially, the background decreases like 
However, the suppression due to ^ and t] in Eq. (3.7) becomes increasingly severe for heavier cj)^. Yet, 



Eq. (p. 61) may have overestimated the signal somewhat due to the omission of a correction factor to 



the capture rate by the sun for a heavier DM particle [12|. The signal rates for the direct and indirect 
searches are shown in Fig. 2 as a function of m^o. The rate for direct search in Fig. 2(a) is for a 
''^Ge detector. The lower solid curve corresponds to a high mass limit for the CP-odd Higgs boson so 
that the dominant contribution is from a SM-like Higgs boson, assuming = 100 GeV. The upper 
solid curve is for an optimal choice of parameters to enhance the signal rate (CP-odd Higgs mass of 
80 GeV and tan/? = 50). In our numerical analysis, the coupling A in Eq. ( p.2| ) has been taken to be 
unity. We see that the rate is typically of order 10~^ events/ (kg day). It is much below the current 
experimental sensitivity of about 2 events/(kg day) [^0|, and is also unreachable by next generation 
experiments with sensitivity of 10^^ events/(kg day) ll^. Fig. 2(b) shows the calculated signal rate 



for the indirect search, with the same sets of SUSY parameters as in (a). Based on Eq. (|3.8D, the 
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background rate is calculated and shown by the dotted curve in Fig. 2(b). Again, it is difficult to 
detect the signal for the indirect search as well. 

4 Conclusions 

We have considered the possibility that the CDM constitutes of sneutrino-like messenger particles. 
We found that the neutral components (j)^ of the SU(2) doublet is naturally lighter than the charged 
component over most of the relevant parameter space due to one-loop electroweak corrections and 
could serve as a CDM candidate. However, our relic density calculation shows that a significant 
amount of messenger CDM in the minimal model has been already ruled out by present experimental 
results. We introduce a mechanism that generates a CP-even-CP-odd mass splitting. This circumvents 
the constraints from present detection experiments and allows the messenger sneutrino to be a valid 
CDM candidate in a wide range of SUSY parameter space. Consequently, the detection via direct and 
indirect searches would be very difficult due to the rather small messenger sneutrino-nucleon scattering 
cross section. 
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Appendix A: m^- m^o mass difference 

At tree-level, the only source of a mass splitting of different members of an SU(2)j;,(8>U(l)y are 
a non-zero vev for the SU(2) D-term which is strongly suppressed as shown in Eq. ( |2.5D . In this 
appendix we present the result for the dominant one-loop contribution to this mass difference. A 
complete one-loop calculation is rather involved. In particular, it requires a renormalization of tan (3 
which enters already at tree-level. Hence, we consider the case tan (3 = 1 where the corrections to the 
mass squared difference is finite and given by the difference of the self-energies 

Am^ = ml- - m|o = yl^-<^-(m|-) - yl^o0o(m^o) • (A.l) 

Furthermore, we neglect terms suppressed by F/M'^ (i. e., we set F = and as a result m^o = m^- = 
M). 

There are three types of diagrams that contribute to the self-energies: the loops involving gauge 
bosons (^f), gauginos {g), and loops with tri-linear terms arising from the D-terms (Higgs field only). 
Loops with quartic terms arising from the D-terms do not generate mass splitting among the different 
members of the multiplet <j). The results are 

(4M2 - ml) Bo{M^,mlM^) - A^im^) - Ao{M^)] , 



a 
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a 



27r sin^ lOw 

—a 



mlBo{M\ml,,M^), 



(A.2) 



47r sin^ 26*11/ 

where a is the fine structure constant and 9w the weak mixing angle. In the SUSY hmit we have 
rriHO = mi = m^. The self-energy for the charged component can be written as 

A.-s-{M^) = cos" 2ewA.oAM^) + sin^ 2ewAso.o{M^){ml ^ 0) 



Hence, we obtain 

Am^ = — [5o(M2, 0, M^) - Bo{M'^,ml, M^) . 
All L 

In dimensional regularization, the scalar one and two point functions are given by 
167r'^Ao{ml) = (a + 1 - Inm^) , 



(A.3) 
(A.4) 



A 



dx In 



xm\ + (1 — x)m\ — x{\ — x)p^ 



(A.5) 



where A = 1/e — 7^+ln47r in d = 4 — 2e dimensions, with 7^ the Euler constant. Expanding Eq. (A.4) 
for M ^ TTiz, Eq. (|2.6| ) is recovered. 



Appendix B: Annihilation Rate and Relic Density 



Our calculation for the annihilation rate follows the formalism developed in Ref. |13|. There are four 
nearly mass-degenerate particles present at the freeze-out temperature Ty, (ff^ and c/)*^^ {Q = 0, — 1). 
If we assume that there is no messenger number asymmetry then the total relic density is n = An^Q = 
4nj,Qt and the rate equation can then be written as 



dn 
di 



-3i7n — (crt'rel)'T- , 



where H is the expansion rate of the universe and 

1 



a = 



Q,Q' 



Note a 



a 



0. 



Following Ref. |13|, we first define a Lorentz invariant function 



wis 



Q,Q,a,b 



ab) 



d cos 9 



(B.l) 



(B.2) 



(B.3) 



where = 1 — 4m^/s is the squared velocity and the sum over a and b symbolizes the sum over all 
two-body final states. Then 



g,g'=b,w 



If 



(B.4) 
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the result of the expansion of the different contributions to w{s) to first order in is 

^4 



w 



9 



(3 - 2/32) , 



w 



w 



bw 



5127r 

(3 - 4/32) , 
g- A 4 



5127r 

■,'4 

5127r V 3''^ 



w 



9 



w 



bw 



W 



W 



..bb 



5127r 
5127r 



12r 



+ 



4/32 



(l+r)2 (l + r)4 
12r 4r/?2 



1 + 7r 
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5127r 



(l+r)2 (l + r)4 
4r 4r/32 / 1 



1 + 6r — 3r' 



(1 +r)2 (1 + r)4 V3 



w 



If 



E 

2w' 



•2 ^ 



Nc0- 
40967r 



iVwK,V' + 6(iVw-l)5' 



(B.5) 



where r = M/m^o. The result for the gauge bosons (gauginos) corresponding U(l)y and S\J{2)l 
are denoted by a superscript 6 and w {b and tt)), respectively. Here, Nc = 1 (3) for colored singlets 
(triplets) and = 1 (2) for SU(2)i singlets (doublets). The index s (/) runs over all MSSM scalar 
(fermion) fields other than gauge bosons and gauginos. 



The thermal average (cj^rci) can be expressed by 

1 



(cr^^rel) 



{A + Bx) , 



where x = Ty/m^o, and |13 | 



A = w{Amlo), B = 3 



(B.6) 



(B.7) 



We have ignored the possible contribution to the annihilation cross sections from the interactions in 
Eq. (|3.2|). This implies a lower limit of {avrei) and hence, an upper limit on the relic density. 



It is customary to express the relic abundance in terms of the mass density in units of the critical 
density = p/pc- It is found that 



8.5 X IQ- 

\f9l 



VTeVy A^\Bx 



(B.8) 



where h is the Hubble constant, ~ 228.75 0]. The value x is related to the freeze-out temperature 
and can be obtained iteratively by 



-1 



m, 



/Tf = In 



0.076 Mpi 



{A + B x) ^/x 



(B.9) 



where Mpi is the Planck mass. 
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Appendix C: sneutrino- nuclei Elastic Scattering via Higgs exchange 



For the model of our interest, the Z-exchange contribution between two different sneutrino mass 
eigenstates is negligible. The leading contribution would therefore be the spin-independent Higgs 
exchange. We write the effective interaction between a sneutrino (0) and a nucleon {N) 



where {N = p, n) is the effective coupling through Higgs boson exchanges: 

u,d,s c,b,t c,b,i 



(C.l) 



(C.2) 



where Kg, Kg and Kq are contributions from a light quark q, a heavy quark Q and its super symmetric 
partner Q. In analogy to the neutralino- nucleon interaction ll^], we determine the couplings as 



m 



m 



2 ' 



108 ■'^'3 



2 2 
3 Q 



(C.3) 



The constant /j.^ is the nucleon mass fraction due to a light quark q, and frp^ = 1 — X^g'^^'^/T,- 
Numerically, we take |21] 



for a proton : 
for a neutron 



/r^= 0.023, /r, = 0.034, 
/r„= 0.019, /7,_^= 0.041, 



(C.4) 
(C.5) 



and 1221 



0.14. 



(C.6) 



The 



couplings are given by 



A V sin j3 cos a, c, 



1 



sin /3 sin a, 



(C.7) 



where A is a model-dependent parameter in Eq. (|3.2| ), v = 246 GeV is the Higgs vacuum expectation 
value and (5 and a are the standard mixing parameters in the SUSY Higgs sector. Finally, g g 
are the couplings of a Higgs boson Hj to quarks and squarks and we have followed the convention in 
Appendix A of Ref. |17]. In practice, the squark contributions are small and have been neglected. 

The cross section for the coherent elastic scattering of i?!> off a nuclei A in the non-relativistic limit 
is thus calculated to be 



m, 



IGvr (m^ -I- niAy 47r (m^ -|- m^)^ 



[fpZ + UiA-Z)f. 



(C.8) 
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